BioData Mining
○ Springer Science and Business Media LLC
Preprints posted in the last 7 days, ranked by how well they match BioData Mining's content profile, based on 22 papers previously published here. The average preprint has a 0.03% match score for this journal, so anything above that is already an above-average fit.
Karalis, V.
Show abstract
The aim of this study was to develop a framework for assessing the incremental value of decisions on drug amount, timing, input structure, cessation and restart. The LOGOS (Longitudinal Optimization of Global Outcomes and Strategies) approach combines matched nested regimen classes, mechanism-specific benchmarks and evidence requirements in a single longitudinal comparison. Thus, the contribution of additional regimen freedom is assessed separately from improvements due to different resources, constraints or evaluation horizons. Nested feasible classes have nondecreasing objective suprema, whereas a strict gain requires a feasible witness above a valid bound for the smaller class. The approach is illustrated through analytical examples of periodic smoothing, matched desensitization, relaxed exposure allocation, finite-course loading and safety after cessation. In a normalized finite Emax model, integrated benefit increased from 0.629606 for three equal boluses to 0.672945 after optimization of amounts, 0.677043 for an amount-and-time candidate and 0.679553 for a globally certified loading-maintenance-coast solution. The adjacent observed gains were 6.88%, 0.61% and 0.37%, respectively, and all schedules satisfied the full-tail injury bound. In addition, a periodic factorial example separated response curvature from a joint recovery and repair process, including a null case and results dependent on feasibility constraints. The LOGOS framework allows assessment of whether additional regimen flexibility improves benefit, provides no added value or restores feasibility. These dimensionless examples support conditional mathematical and model-based conclusions; clinical superiority and safe restart in patients have not been established.
Nguyen, D. H.; Pham, H. T. V.; Ding, J.; Tarning, J.
Show abstract
Aims Rosuvastatin is commonly co-administered with ticagrelor in patients with acute coronary syndrome. Ticagrelor inhibits breast cancer resistance protein (BCRP/ABCG2), increasing rosuvastatin exposure and potentially the risk of rhabdomyolysis. This risk may be further influenced by age and the ABCG2 c.421C>A variant. We aimed to predict optimal rosuvastatin dosing with ticagrelor in elderly patients stratified by ABCG2 genotype using physiologically based pharmacokinetic (PBPK) modelling Methods PBPK models for rosuvastatin and ticagrelor were developed and validated in PK-Sim using clinical data from healthy adults, elderly patients, and individuals with wild-type or homozygous ABCG2 c.421C>A genotypes. Key parameters, including BCRP catalytic activity and age-dependent declines in CYP2C9 and CYP3A4 activity, were optimized. The drug-drug interaction model incorporated competitive BCRP inhibition by ticagrelor. Rosuvastatin steady-state exposure was simulated with half or quarter standard doses during ticagrelor co-administration, and exposure matching approach was used to identify optimal dosing. Results The PBPK models demonstrated good predictive performance across populations and genotypes. In young adults with the wild-type genotype, halving the rosuvastatin dose resulted in an AUC fold change of 1.44, supporting a half-dose strategy. In young homozygous variant carriers, half-dose co-administration produced comparable exposure (AUC fold change 0.78). In elderly wild-type individuals, a four-fold dose reduction achieved comparable exposure (AUC fold change 1.06), whereas elderly variant carriers required only a half-dose reduction (1.12). Conclusions PBPK modelling supported age- and ABCG2 genotype-specific rosuvastatin dosing with ticagrelor. Elderly wild-type patients may require a further dose reduction to minimise toxicity while maintaining therapeutic exposure.
Lambert, B. L.
Show abstract
Background. Regulatory pressure to avoid confusingly similar drug names appears to have pushed manufacturers to coin more distinctive proprietary names, but no one has quantitatively analyzed this trend. Objective. To determine whether proprietary drug names became more orthographically distinctive (unusual in spelling) over time. Methods. We examined 1,321 US Food and Drug Administration (FDA) new molecular entity (NME) names (1985-2024), 2,769 FDA Orange Book trade names (1982-2026), 177,337 US pharmaceutical trademarks (1980-2026), and 1,273 generic names (1982-2025). Compared to prior names (including those already marketed) and ordinary English, we scored each name's character trigram (three-letter sequence) surprisal (statistical unexpectedness) and its edit-distance similarity to the closest existing names. Linear regression tested time trends, controlling for name length. Results. Mean NME trigram surprisal rose in every decade between the 1980s and 2020s, from 7.21 to 8.50 nats (natural-logarithm units, higher meaning less predictable), an 18% increase (b = 0.038, p < .001, R2 = 0.34). Half of 2020s names (7% in the 1980s) contained a letter sequence found in no common English word (e.g., Jesduvroq, Omvoh). Nearest-neighbor similarity fell after the early 1990s. Trademark data showed the same trend (b = 0.016, p < .001). Generic names showed a weaker trend (b = 0.011, R2 = 0.063). Conclusions. Over four decades, drug names became more orthographically distinctive, consistent with namespace crowding, deliberate efforts to avoid confusing names, and branding strategy. Because unfamiliar letter sequences can introduce other error pathways, increased distinctiveness does not necessarily imply improved safety.
Workalemahu, T.; Nicholas, T. J.; Quinlan, A. R.; Yang, X.; Griffin, L. B.; Leisher, S.; Blue, N. R.; Arthurs, A. L.; Page, J. M.; Branch, D. W.; Camp, N. J.; Jorde, L. B.; Roberts, C. T.; Jessica, G. L.; Wapner, R.; Kingdom, J.; Gibbins, K. J.; Parast, M. M.; Reddy, U. M.; Silver, R. M.
Show abstract
Abstract Background: Placental insufficiency (PI), characterized by impaired placental function, is a leading cause of stillbirth. Although rare coding variants are known to cause a subset of stillbirths, their contribution to PI is poorly understood. Methods: We conducted a secondary analysis of stillbirth cases (N=292) of the Stillbirth Collaborative Research Network study with exome sequencing available on postmortem placental or fetal samples. PI was defined by standardized placental pathological lesions (e.g., maternal vascular malperfusion) and prespecified clinical criteria for placental dysfunction. Candidate gene-based rare-variant burden analysis was performed using the Sequence Kernel Association Test-Optimal framework, focusing on evolutionarily constrained developmental and placental-function genes. Primary analysis compared stillbirths with and without placental vascular malperfusion. Similar comparisons by clinical PI and other placental pathological lesions were conducted in secondary analysis. Results: Stillbirths with PI phenotypes demonstrated different patterns of genetic variation than stillbirths without PI. Specifically, we observed nominal differential enrichment of constrained variants in biologically prioritized placental-function genes, regulating extravillous trophoblast invasion and spiral artery remodeling (CTNND2, PHACTR1, SLIT1, VEGFA), placental angiogenesis (IGF2R, TEK), mitochondrial homeostasis (OPA1), and fetal growth (DLK1, GNAS). Secondary analyses on specific placental pathological lesions suggested recurrent enrichment of genes within similar placental developmental pathways. Conclusions: We observed nominal differential enrichment of constrained developmental and placental-function genes in stillbirths with PI compared to stillbirths without PI. The observed genetic signals across PI phenotypes suggest shared biological pathways underlie placental dysfunction. Our findings advance genotype-phenotype characterization and help prioritize candidate pathways for larger collaborative genomic and functional studies to improve stillbirth risk stratification and prevention.
Balakrishna, K.; Bhatt, V.; Hammond, A.; Bawoke, M. M.; Fradley, M. G.; Zlotoff, D. A.; Neilan, T. G.
Show abstract
Background. Immune checkpoint inhibitors (ICIs) can cause immune-related adverse events in almost any organ. While myocarditis has received the most attention, data on pericardial disease with ICI use are limited. We examined pericardial adverse event reporting with nine ICIs using the FDA Adverse Event Reporting System (FAERS). Methods. We combined legacy AERS and current FAERS quarterly files from January 2011 to June 2026 into 18,629,240 unique cases, deduplicated by case number, and removed cases the FDA had listed as deleted. ICI exposure was identified from generic names, brand names, development code names and the active-ingredient field. Pericardial events were defined by seven MedDRA preferred terms. We estimated reporting odds ratios (ROR), proportional reporting ratios (PRR), Bayesian information components (IC) and empirical Bayes geometric means (EBGM). We also examined myocarditis co-reporting, time trends, a comparison with other antineoplastic agents, a MedDRA query-based outcome, and 14 alternative definitions. Results. The analysis included 18,514,045 cases: 277,434 involved an ICI, 34,792 carried a pericardial term, and 1,680 had both. Pericardial events were reported disproportionately with ICIs (ROR 3.35 (95% CI 3.19 to 3.52); PRR 3.34; IC025 1.62). Eight of nine agents had a ROR above 1 (3.32 to 5.63). Among the most frequently reported terms, estimates were highest for cardiac tamponade (5.59 (95% CI 4.98 to 6.28)), followed by pericardial effusion (3.83 (95% CI 3.60 to 4.07)) and pericarditis (2.09 (95% CI 1.89 to 2.31)). After removing cases that also reported myocarditis, the ROR was 3.13 (95% CI 2.98 to 3.30). Most sensitivity analyses gave similar estimates (ROR 3.13 to 3.47); defining exposure by generic drug names alone, which identified only 44% of exposed cases, gave 3.99. Compared with cases naming other L01 antineoplastic agents, the ROR was 1.13 (95% CI 1.08 to 1.19): 1.28 (95% CI 1.20 to 1.36) for ICI-only cases and 0.96 (95% CI 0.89 to 1.04) when another antineoplastic was also named. The ICI-specific reporting odds ratio fell each quarter (odds ratio 0.977 per quarter, 95% CI 0.973 to 0.980). Conclusions. In FAERS, pericardial events were reported more often with ICIs than with other drugs across most agents and terms; however, incidence and causation cannot be determined. The excess over other antineoplastic agents was small, and reporting disproportionality declined over time. These findings support clinical awareness of pericardial disease in patients receiving ICIs and warrant study in data that can measure incidence.
Olivier, L.; Madzhakova, D.; Lauterbach Oprea, C.; Stollenwerk, A.; Orlikowsky, T.; Schoberer, M.
Show abstract
Introduction: Fighting the ventilator (FV) is characterized by a patients active counteraction of respiratory gas flow across multiple breaths during mechanical ventilation. An effective ventilation is prevented during such episodes. The objective of this exploratory study was to evaluate clinicians knowledge, perception and agreement regarding this condition. Methods: In a cross-sectional online survey, medical professionals with experience in intensive care medicine (mostly neonatal and pediatric patients) responded to 11 questions covering various clinical aspects of FV. Results: All respondents were familiar with the phenomenon and assessed it to be clinically relevant. Neonates were perceived to be predominantly affected. Promoting factors of FV included inadequate sedation and ventilation-associated factors. In clinical practice, the participants identified FV by clinical evaluation, analysis of ventilator waveforms and capnography and used both pharmacological and non-pharmacological intervention strategies. There was no consensus concerning the English terminology of the condition. Only few respondents were familiar with local guidelines to manage FV. Conclusion: Although there was consensus on the relevance of FV, diagnostics and management strategies were highly heterogeneous. Many participants agreed that poor patient-ventilator synchrony was an important cause of FV. Future work should focus on characterizing and generating evidence on FV.
Lefebvre, A.; Benusiglio, P.; Nuel, G.
Show abstract
Background: Current Lynch syndrome (LS) prediction models do not account for detailed family structure and tumour mismatch repair (MMR) status. We developed LynchRisk, a Bayesian model that estimates individualized LS probability from pedigree, cancer history and tumour MMR test results. Methods: LynchRisk combines Mendelian inheritance, age- and sex-specific cancer hazards, affected and unaffected relatives, and detailed MMR status. Parameters were derived from the literature. Results: In two simulated cases of probands with LS-associated cancers, family structure and MMR status markedly changed LS probability estimates, with implications for germline testing decisions. LynchRisk also identified situations in which a negative germline result might be considered sufficient to rule out LS or, conversely, prompt further investigations such as somatic MMR or long-read sequencing. Discussion: LynchRisk gives individualized LS probability by combining clinical, family, and tumour data at different stages of the diagnostic process. Further work will focus on recalibrating cancer hazard functions, validating model performance through simulations, and developing an R package for use by interested research groups.
Kagawa, E.; Kato, M.; Hidaka, T.; Kunita, E.; Hironobe, N.; Matsui, S.; Kobayashi, Y.; Uotani, Y.; Teramoto, T.; Yoshii, K.; Okumoto, T.; Tani, M.; Dote, K.
Show abstract
Background: Artificial intelligence (AI) may improve electrocardiographic differentiation of Takotsubo syndrome from acute anterior ST-segment elevation myocardial infarction (STEMI), but its impact on human diagnostic performance, the added value of visual explanations, and differences by reader background and experience remain uncertain. Methods: Using 12-lead electrocardiogram (ECG) images from two hospitals, we developed and externally validated an ensemble of three deep learning models after self-supervised pretraining and fine-tuning. The diagnostic task used ECGs recorded within 14 days after symptom onset, regardless of reperfusion status. We conducted a four-phase reader study. Phases 1-3 compared unaided interpretation, a reference sheet, and additional AI assistance. In Phase 4, readers evaluated a new ECG set after randomization to reference-sheet support alone, with AI predictions, or with AI predictions plus gradient-weighted class activation mapping (Grad-CAM) heatmaps. Diagnostic correctness was analyzed using generalized linear mixed-effects models. Exploratory analyses examined differences by cardiologist status and experience. Results: Development and external-validation cohorts, including healthy controls, comprised 832 and 95 participants, respectively. For Takotsubo-STEMI discrimination, external-validation accuracy was 0.914, with areas under the receiver operating characteristic and precision-recall curves of 0.970 and 0.984. Seventy healthcare providers completed the reader study. The estimated probability of correct diagnosis increased from 0.632 without assistance to 0.864 with AI assistance (odds ratio [OR], 3.71; P<0.001). In randomized Phase 4, AI assistance improved accuracy vs reference-sheet support alone (OR, 4.28; 95% CI, 3.08-5.95; P<0.001), and Grad-CAM provided additional benefit (OR, 1.53; 95% CI, 1.04-2.25; P=0.033). Exploratory analyses showed greater improvement among less experienced readers; with AI assistance, non-cardiologists outperformed cardiologists. AI-assisted readers remained less accurate overall than the model alone. Conclusions: An ECG image-based ensemble differentiated Takotsubo syndrome from anterior STEMI and improved reader accuracy. Visual explanations provided additional benefit, while experience-related differences highlight the importance of effective human-AI interaction.
BOU ALI, M. N.; Mahmoud, M.
Show abstract
Background. Real-world 12-lead electrocardiograms (ECGs) often have missing or corrupted leads. A deep-learning classifier can be prepared for this by random lead masking during training, or by injecting clinical knowledge, for example by pulling its attention towards the leads that guidelines use for each diagnosis. Whether this knowledge adds robustness beyond masking is unknown. Methods. Pre-registered comparison of three arms sharing one network (per-lead 1D-ResNet encoder and per-diagnosis attention over the 12 leads): standard training (M0); random whole-lead masking (M1); masking plus a Kullback-Leibler penalty pulling attention towards a fixed guideline-derived key-lead map (M2, guideline-based attention guidance). Six diagnoses; five seeds per arm. Development on PTB-XL (21,799 ECGs); single sealed evaluation on PTB-XL fold 10 (2,198 ECGs) and on CODE-15% (233,204 exams, external validation). Primary endpoint: AUROC difference M2 - M1, averaged over first-degree AV block, right and left bundle branch block and atrial fibrillation and over 1 to 6 lost (zeroed) leads, with a patient-level paired BCa bootstrap 95% confidence interval. Results. The guidance did not change robustness: {Delta}AUROC -0.0010 (95% CI -0.0033 to +0.0006), an interval entirely within the pre-specified +-0.02 margin; the result was unchanged with another bootstrap seed, other interval methods and leave-one-diagnosis-out analyses. Masking improved AUROC: M1 - M0 = +0.0134 (95% CI +0.0102 to +0.0172) when one key lead was lost, and the number of lost leads tolerated before AUROC fell by more than 0.02 rose from 7 (M0) to 10 (M1). M2 attention matched the map and M2 was non-inferior on intact ECGs, but its calibration slope was lower for all six diagnoses in both datasets. M2 lost slightly more when all key leads were lost together (-0.0067) and gained only when one to three leads remained (exploratory; AUROC 0.9408 vs 0.9210 with 11 leads lost). CODE-15% reproduced every pattern (primary -0.0007, 95% CI -0.0011 to -0.0003). Conclusions. Random lead masking suffices to make an ECG classifier robust to missing leads; adding guideline knowledge aligns attention with clinical criteria but adds no robustness and worsens calibration. Guideline-aligned attention should not be taken as evidence of robustness.
Yamada, E.; Stark, A.; Burdine, J.; Donikini, L.; Datta, P.; Krutsch, K.
Show abstract
Objective: Amphetamine (AMP) and its derivates are commonly prescribed for the management of attention-deficit/hyperactivity disorder (ADHD). This study aimed to 1) quantify AMP's transfer into human milk in lactating women receiving therapeutic doses and 2) contextualize infant exposure. Methods: Milk samples and dyad health histories of thirteen lactating women treated with varying doses of either mixed AMP, dextroamphetamine, or lisdexamfetamine were released for this study from the InfantRisk Human Milk Biorepository. Milk drug concentrations were quantified by validated liquid chromatography-tandem mass spectrometry (LC-MS/MS) and used to calculate standard relative infant dose (RIDmat) and RID standardized to the defined daily dose (RIDDDD) of 15 mg AMP base. Results: Average milk concentrations (Cavg) of AMP base ranged from 23.3 to 297.0 ng/mL and increased linearly (R^2 = 0.92) with maternal daily dose. The RIDmat ranged 3.6-10.8% and RIDDDD ranged 1.4-21.8%. The fitted regression line crossed the 10% RIDDDD threshold at 1.17x DDD, with the upper 95% prediction limit crossing at 0.84x DDD. No adverse event reported by the mother-infant dyads were attributed to AMP exposure, notably including absence of reduced milk supply. Conclusions:The levels of AMP observed in human milk followed a dose predictable manner. The average patient was predicted to remain below the 10% RID DDD safety threshold at total daily doses at or below approximately 28 mg of mixed AMP salts or 59 mg of lisdexamfetamine. Ninety-five percent of patients were predicted to remain below the threshold at or below approximately 20 mg and 40 mg, respectively.
Baierl, J. D.; Tyrer, J. P.; Barnes, D. R.; Jones, M. R.; Peng, P.-C.; Tischkowitz, M.; Antoniou, A. C.; Pharoah, P. D. P.
Show abstract
Polygenic scores (PGSs) have progressed toward clinical utility in recent years, demonstrating odds ratios that are comparable with both lifestyle-related and monogenic risk factors for several complex diseases. However, PGS implementation relies on developing analytically and clinically valid assays. Methods for developing polygenic models (PGMs) that incorporate tens of thousands to millions of variants require either genome-wide chip genotyping combined with imputation or whole genome sequencing, limiting their portability. Additionally, variants may be in regions where the local genomic sequence makes it difficult to assay, and multi-allelic or strand ambiguous variants can make bioinformatics more complicated. Because of these challenges, attempts to develop a reliable assay for published PGMs can produce multiple variations from the same underlying model, hindering routine clinical deployment and incorporation into multi-factorial risk prediction models. In this analysis, we present an adjustment to the previously published select and shrink with summary statistics (S4) PGM method designed to optimize performance for realistic clinical use by restricting to a pre-specified number of available SNPs that are well-captured by genotyping arrays. We compare the performance of this oligogenic S4 method (S4-Oligo) to currently published PGSs as well as the full S4 model for prediction of breast, colorectal, ovarian, and prostate cancers in the All of Us Research Program, Biobank Japan, and UK Biobank. We find that S4-Oligo matches or exceeds the accuracy of published PGMs at comparable model sizes while retaining between 88% and 99% of the performance of the original S4 model when limited to 1,000 variants.
Celik, M. C.; Cicek, Y.; Sakiroglu, M. A.; Aslan, D.; Altunc, A. T.; Velioglu, H. A.
Show abstract
The brain differences of dyslexia are well described and its genetic risk is now measurable; the brain systems that connect the two are unknown. That link decides which measures complement one another in early identification and which systems intervention research should examine. Here we compared large-scale polygenic score maps of brain structure in dyslexia with normative maps of brain neurotransmitter systems. The structural signature of genetic liability aligned with the cortical distribution of neurotransmitter systems. The strongest and most consistent alignment appeared for the histamine H3 receptor, the vesicular acetylcholine transporter (VAChT) and the N-methyl-D-aspartate (NMDA) glutamate receptor. The alignment was most pronounced in cingulate, medial and superior frontal cortex and the paracentral lobule. The H3 alignment held after adjustment for cortical position and for metabolic and cognitive reference maps. Once the component common to all nineteen maps and position were removed together, no individual system remained significant, while the shared component itself aligned with liability. The signature therefore follows the topography these systems share. In 169 children, per-person alignment with the H3 map related to reading skill with the opposite sign to liability. To our knowledge, this is the first study to relate the structural signature of genetic liability for dyslexia to normative maps of neurotransmitter receptor and transporter density. It places that signature within the shared geography of the systems that set attention, arousal and synaptic plasticity during learning, and gives the genetics of dyslexia a testable biological axis.
Khan, J. Y.; Groh, J. M.; Dookeran, K. A.
Show abstract
Objectives: Using a nationally representative database of U.S. hospitalizations, we aimed to (1) characterize contemporary trends in hypoxic-ischemic encephalopathy (HIE) prevalence by severity, therapeutic hypothermia (TH) utilization, and associated in-hospital mortality; and (2) examine the association between TH and mortality by HIE severity and gestational age (GA), focusing on 35- and 36-weeks GA. Methods: We analyzed newborn hospitalizations from the National Inpatient Sample, 2016-2023. Unspecified, mild, moderate, and severe HIE were identified using ICD-10-CM diagnosis codes P91.60-P91.63; and 35- and 36-weeks GA using codes P07.38-P07.39. TH was identified using ICD-10-PCS codes 6A4Z0ZZ and 6A4Z1ZZ. Survey-weighted multivariable generalized linear models with marginal standardization were used to estimate adjusted mortality risk differences (RD) with surprisal (S) values, and assess whether associations between TH and mortality differed by HIE severity across GA strata. Results: Among 29,767,238 hospital discharges, HIE was identified in 66,535 neonates (2.24 per 1,000), and TH used in 19.2% of HIE cases. In-hospital mortality was 9.74% among all HIE cases and 6.65% among TH-treated neonates. From 2016-23, HIE prevalence nearly doubled (1.58 to 3.01 per 1,000; p-for-trend <0.001) and TH utilization more than doubled (0.33 to 0.74 per 1,000; p-for-trend <0.001), while TH use as a proportion of HIE remained stable (17.8%-20.8%; p-for-trend=0.39). HIE severity distribution was: mild (15%), moderate (29.4%), severe (18.3%), and unspecified (37.3%), with moderate HIE increasing over time (22.3%-33.5%). Crude mortality was consistently highest for severe HIE (range: 33.65%-35.25%). Among infants with severe HIE, TH was associated with reduced mortality overall (RD -10.12%; 95% CI -14.86, -5.39; p<0.001; S=15.12), and specifically at 36-weeks GA (RD -21.98%; 95% CI -35.85, -8.12; p=0.002; S=9.06), but not at 35-weeks (RD -1.83%; 95% CI -21.2, 17.55; p=0.86; S=0.22). Conclusions: TH was associated with mortality reduction in severe HIE overall and at 36-weeks GA, but this benefit remains uncertain at 35-weeks GA.
Sun, Y.; Morgante, F.; Mackay, T.
Show abstract
Polygenic scores (PGS) predict an individual's genetic value for a complex trait of interest by aggregating effect size estimates from genome-wide association studies (GWAS). While powerful, this approach assumes that genetic variants affect the phenotype in an additive way, which might be problematic when dominance is part of the genetic architecture. To address this limitation, we developed PGSa,d, a simple polygenic scoring framework, rooted in quantitative genetics theory, that incorporates both additive (a) and dominance effects (d) directly. We evaluated PGSa,d against the standard GWAS-based PGS (PGS{beta}) in simulations. The results showed that PGSa,d consistently achieved lower root mean squared error (RMSE) and progressively higher R2 as dominance contributed more to the genetic architecture. We then applied the method to 16 blood cell traits, height, and body mass index (BMI) from UK Biobank. PGSa,d showed comparable R2 to PGS{beta} for all traits (including height and BMI whose architecture is predominantly additive) and significantly lower RMSE for four traits after strict Bonferroni correction. Across traits, RMSE improvement from PGSa,d was positively associated with the magnitude of the estimated dominance effects, whereby traits with larger dominance effects displayed lower RMSE from PGSa,d. As a byproduct of our analysis, we also identified 141 variants in 39 genomic regions with significant dominance effects after Bonferroni correction across traits, recapitulating known associations and identifying novel ones. Together, these results show that explicitly modeling additive and dominance effects is a robust strategy that can improve individual-level predictions, especially when dominance contributes to complex trait architecture.
Sharew, N. T.; Clark, S. R.; Arolt, V.; Brockmoller, J.; Chen, C.-H.; Domschke, K.; Hall-Flavin, D. K.; Illi, A.; Kampman, O.; Leinonen, E.; Liou, Y.-J.; Mushiroda, T.; Stingl, J. C.; Benjaponpitak, A.; Mahasirimongkol, S.; Wang, L.; Kato, M.; Liu, Y.-L.; Bobo, W. V.; Tsai, S.-J.; Kubo, M.; Weinshilboum, R. M.; Mrazek, D. A.; Biernacka, J. M.; Baune, B. T.; Schubert, K. O.; Amare, A. T.
Show abstract
Biologically informed polygenic scores (PGSs) may help interpret genetic associations with treatment outcomes. Here, we computed neurotransmitter-based PGSs (NeuroPGSs) for patients with major depressive disorder (MDD) and assessed their association with SSRI treatment outcomes. By targeting genetic variants involved in key neurotransmission pathways, we developed PGSs for serotonergic (5228), gamma-aminobutyric acid (GABA, 4103), dopaminergic (7061), and glutamatergic (5690) variants. We leveraged data from two cohorts of patients who received SSRIs for at least 4 weeks: the International SSRI Pharmacogenomics Consortium (ISPC) and the Pharmacogenomics Research Network Antidepressants Medication Pharmacogenomics Study (PGRN-AMPS). Treatment response was defined as [≥]50% reduction from baseline depression symptoms, measured using 17-item Hamilton Rating Scale for Depression (HRSD-17) in ISPC and 16-item Quick Inventory of Depressive Symptomatology (QIDS-C16) in PGRN-AMPS. Remission was defined as HRSD-17 [≤]7 (ISPC) or QIDS-C16[≤]5 (PGRN-AMPS). The association between each NeuroPGS and antidepressant response or remission was evaluated using logistic regression models, adjusting for age, sex and five genetic principal components, with statistical significance determined at Benjamini Hochberg-adjusted p<0.05. The GABAergic-based PGS was associated with SSRI response (aOR per SD = 1.45, 95%CI: 1.14-1.85). Patients with the highest GABAergic PGS (10th decile) had 2.59-fold higher odds of response (95%CI: 1.46-4.60) compared to those with the lowest decile (1st decile). These exploratory findings support further evaluation of NeuroPGSs; independent replication and assessment of incremental predictive value are needed before clinical application. Keywords: Polygenic score; Neurotransmitters; GABAergic; depression; SSRIs; pharmacogenomics
Hew, C. K. W.; Leserri, S.; Aydogan, D. B.; Bastin, M. E.; Bonthrone, A. F.; Paterson, K.; Barclay, G.; Jardine, C.; Smith, K. M.; Thrippleton, M. J.; Marioni, R. P.; Boardman, J. P.; Blesa Cabez, M.
Show abstract
Structural connectomics provides a means of characterising white matter (WM) network organisation using diffusion MRI (dMRI), but reconstructing neonatal structural connectomes remains challenging due to the immature and heterogeneous nature of the developing brain. Although recent neonatal studies have increasingly adopted approaches beyond conventional diffusion tensor imaging (DTI), there remains scope to incorporate more advanced tractography and streamline post-processing approaches within neonatal pipelines. Here, we propose a neonatal structural connectome pipeline integrating surface-based anatomical constraints, ensemble tractography, hierarchical anatomical refinement of pathways (HARP) filtering, and spherical-deconvolution informed filtering (SIFT2) reweighting. We acquired 3T MRI data from 300 neonates (185 preterm, 115 term-born) to assess the impact of pipeline configuration on network organisation. Mean gestational age (GA) at birth was 29.31 (23.43-32.86) weeks in preterm infants and 39.52 (37.00-42.43) weeks in term infants. At scan, mean GA was 40.85 (36.29-49.00) weeks in preterm infants and 42.06 (38.29-48.00) weeks in term infants. Four connectome variants were generated: raw surfaced-based tractograms, HARP-filtered, SIFT2-weighted, and combined HARP+SIFT2. Pipeline performance was primarily evaluated using an existing brain age prediction model: connectome-based predictive model (CPM) and our derived connectomes to estimate GA at scan in term and preterm infants. In addition, global graph theoretical metric analyses were conducted to provide a complementary assessment of pipeline configuration differences in network organisation. In total, seven graph theoretical metrics were calculated and evaluated: global efficiency, local efficiency, transitivity, modularity, assortativity, degree variance and statistical complexity. Across all configurations, connectomes supported accurate brain age prediction, with mean absolute errors (MAE) of 0.76-0.79 weeks in term infants and 1.00-1.06 weeks in preterm infants. Prediction error and brain age gap variability were greater in preterm than term infants across configurations, with HARP+SIFT2 achieving the lowest preterm MAE (1.00 weeks) and brain age gap (BAG) variability (SD = 1.38 weeks). However, differences in prediction performance across configurations were minimal, suggesting limited influence of pipeline choice on brain age prediction. However, analysis of global network metrics demonstrated that pipeline choice significantly influenced downstream measures of network organisation, with HARP+SIFT2 being the only configuration to identify significant preterm vs. term differences across all seven metrics examined. These findings demonstrate the importance of considering tractography and streamline-processing choices when constructing neonatal structural connectomes and provide a framework for more anatomically informed investigation of early brain development.
Noh, S. A.; Park, J.; Yoon, Y. E.; Jeon, J.; Jang, Y.; Lee, S.-A.; Bak, M.; Choi, H.-M.; Hwang, I.-C.; Cho, G.-Y.
Show abstract
Aims: Myocardial strain enables sensitive assessment of cardiac function, but its clinical application has largely focused on individual chambers or selected patient populations. We evaluated the feasibility and clinical relevance of fully automated multi-chamber strain assessment in consecutive patients undergoing routine transthoracic echocardiography (TTE). Methods and results: This retrospective study included 8,350 consecutive adults undergoing TTE. A hospital-wide artificial intelligence (AI)-based system automatically measured left ventricular (LV) global longitudinal strain (LVGLS), left atrial (LA) reservoir strain (LASr), and right ventricular (RV) free-wall strain (RVFWS). AI-derived LVGLS was available from at least one apical view in 97.1% of patients and from all three views in 94.8%; LASr and RVFWS were available in 94.1% and 94.5%, respectively, resulting in complete multi-chamber strain assessment in 90.3%. Strain measures showed modest correlations with conventional chamber-specific echocardiographic parameters and differential associations with clinical outcomes. Gaussian mixture modeling identified three distinct clusters characterized by preserved multi-chamber function, RV-dominant dysfunction, and pan-chamber dysfunction. Compared with the preserved cluster, pan-chamber dysfunction was associated with the highest risk of all-cause death or heart failure hospitalization (adjusted hazard ratio, 3.44; 95% CI, 2.43-4.88; P<0.001). Prognostic associations remained evident in patients with LV ejection fraction (LVEF) [≥]50%, and cluster assignment provided incremental prognostic information beyond clinical factors, LVEF, and RV fractional area change. Conclusions: Fully automated multi-chamber strain assessment was feasible in a large consecutive TTE cohort. Integrating strain across cardiac chambers may provide complementary functional characterization and identify clinically relevant clusters with prognostic information beyond conventional assessment.
Mason, K. J.; Jordan, K. P.; Thomas, M. J.; Bailey, J.; Mughal, F.; Protheroe, J.; Jowett, S.; Saxne Jöud, A.; Snell, K. I.; Dunn, K. M.
Show abstract
Objective: to estimate reconsultation risk and trajectories after a new episode of musculoskeletal pain in children and young people (CYP), and to determine routinely recorded prognostic factors associated with reconsultation. Design: Population based retrospective cohort study. Setting: Clinical Practice Research Datalink (CPRD) Aurum, 2005-2021. Participants: CYP aged 8-18 years consulting primary care for a new musculoskeletal pain episode (no musculoskeletal consultation in prior 12-months). Main Outcome Measures: Musculoskeletal pain reconsultation within 12-months, and 3-year musculoskeletal reconsultation trajectories derived using latent class growth analysis. Prognostic factors selected from a literature search and two-year symptom histories (sample prevalence >1%) were examined using regression models (conditional logistic; multinomial). Results: Among 313,830 CYP with >12 months' follow-up, 32% reconsulted for musculoskeletal pain within 12 months. Among 240,384 CYP with >3 years' follow-up, three trajectories were observed: single episode (62%), early reconsultation (18%), and persistent reconsultation (20%). Older age; female gender; hip, knee, multisite and back pain on index; previous headache; previous abdominal pain; and analgesic prescribing were associated with higher reconsultation risk. Asian or Black ethnicity, and residing in more deprived areas, were associated with lower recorded reconsultation. Conclusion: Most CYP had no further recorded consultation, but 38% followed early or persistent reconsultation trajectories. Routinely recorded prognostic factors, particularly index pain side and prior pain/symptom history, may differentiate groups who could benefit from proactive follow-up. Lower recorded reconsultation among Asian and Black individuals, and living in more deprived areas, should not imply lower need and warrants equity-sensitive monitoring of access, health-seeking, and coding.
Allsopp, S.; Steare, T.; Rojas Arangoitia, V.; Milkias, B.; Ellanki, R.; Dang, H.-A.; Demkowicz, O.; Patalay, P.
Show abstract
Background: Girls are more likely than boys to experience psychological distress during adolescence. Drivers of this gender gap remain under-studied, despite its large impact on later life and outcomes, especially for women. Explanatory factors can be gendered in their prevalence (differential exposure) and/or in their impact (differential susceptibility). We developed a multi-domain framework of explanatory factors for distress and applied it in four low- and middle-income countries. Methods: We first proposed a framework of explanatory factors for psychological distress grouped into seven domains (biological, cognitive, psychological, behavioural, social, environmental and broader societal). We applied this framework to identify explanatory variables assessed in childhood and adolescence in the Young Lives longitudinal cohort studies in Ethiopia (N=908), India (N=952), Peru (N=634) and Vietnam (N=887). Distress symptoms at 18-19 years were measured with the Strengths and Difficulties Questionnaire Emotional Symptoms scale. Gendered prevalence of each explanatory factor was tested with Mann-Whitney tests, and gendered impact with country-specific linear regressions including a gender-by-factor interaction term. Exploratory boosted-tree machine-learning models, interpreted with SHAP values and permutation feature importance, examined the contribution of these factors to the gender gap in each country. Results: Distress was higher in girls than boys in every country, the gap ranging from 0.62 points in Ethiopia and Vietnam to 1.76 in India and 2.45 in Peru. Some explanatory factors were consistently gendered in prevalence (e.g. girls experienced greater role overload in each country), while others were gendered only in some countries (e.g. study amount was higher for girls in Peru and Vietnam, but not Ethiopia or India). Gendered impacts on distress displayed more variability across countries. For example, higher self-esteem was more protective for boys than girls in India, Peru and Vietnam, but more protective for girls in Ethiopia. The machine-learning models identified important predictors across all framework domains, though their salience varied by country. Conclusions: Our findings underscore the importance of considering multi-domain explanations for mental-health outcomes and observed inequalities in mental health. Large cross-country differences in the gender gap and in the prevalence and influence of explanatory factors caution against generalising drivers of mental-health outcomes and inequalities across contexts.
Schuck, R. K.; Moura, I.; Baker-Ericzen, M. J.; Yang, L.-Q.; McDonald, K.; Kapp, S.; Uljarevic, M.; Love, J.; Vera, J. W.; Horner-Johnson, W.; Kripke-Ludwig, R.; Siddeek, Z.; Nicolaidis, C.
Show abstract
Autistic adults face considerable challenges in receiving adequate support, from both paid/formal (such as from disability services) and unpaid/informal (for example, from friends and family) sources. Having adequate support can lead to better health and well-being outcomes. However, most instruments designed to assess adequacy of support were not made specifically for/with autistic adults. This paper presents a psychometric validation of the AASPIRE Adequacy of Supports and Information Scale (AASIS), a community-informed adaptation of a previously developed instrument (the Quality of Care and Support scale; Lucas-Carrasco et al., 2011), designed to be relevant and accessible to autistic adults. We also used the AASIS to explore correlates of perceptions of support in a sample of 798 autistic adults. Our confirmatory factor analysis supported a two factor structure for the AASIS (factor one: Supports; factor two: Information). Further testing confirmed the instrument was invariant across several groups, including race/ethnicity, gender, intellectual disability, communication fluency, and self versus proxy-report. Regression analyses suggested that the way in which participants joined the study (i.e., through community advertisements or via established healthcare/disability services) impacted whether they viewed their support as adequate. Additionally, a higher number of adverse childhood events and greater experience of discrimination were related to worse perceptions of support adequacy across both factors. Our findings support the preliminary validity of the AASIS across a variety of autistic people, though further testing is warranted. Findings also point to the importance of addressing the needs of autistic people who are minoritized and/or who have complex trauma histories.